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Abstract 

The spin symmetry in the Dirac negative energy spectrum and its origin are investigated for 
the first time within the density-dependent relativistic Hartree-Fock (DDRHF) theory. Taking 
the nucleus 16 as an example, the spin symmetry in the negative energy spectrum is found to 
be a good approximation and the dominant components of the Dirac wave functions for the spin 
doublets are nearly identical. In comparison with the relativistic Hartree approximation where 
the origin of spin symmetry lies in the equality of the scalar and vector potentials, in DDRHF 
the cancellation between the Hartree and Fock terms is responsible for the better spin symmetry 
properties and determines the subtle spin-orbit splitting. These conclusions hold even in the case 
when significant deviations from the G-parity values of the meson-antinucleon couplings occur. 

PACS numbers: 21.10.Hw, 21.10.Pc, 21.60.Jz, 24.10.Jv 
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I. INTRODUCTION 



The relativistic Hartree approximation or relativistic mean field (RMF) theory jjj has 
received much attention due to its successful description of infinite nuclear matter as well 
as finite nuclei near and far away from the (3 stability line [24]. One of its great success 
is the natural description of the nuclear spin-orbit potential, which leads to a remarkable 
spin-orbit splitting for the states with the same orbital angular momentum and opposite 
spins (J = I ± 1/2), allowing for the understanding of the magic numbers and forming the 
basis of nuclear shell structure. Furthermore, the pseudo-spin symmetry [5, (|, i.e., the 
near degeneracy between two single-particle states with the quantum numbers (n,l,j = 
I + 1/2) and (n — 1, 1 + 2, j — I + 3/2), whose origin was a long mystery in nuclear physics 
[3, Isj], is well interpreted within the relativistic scheme with local potentials (see Ref. 9] 
and references therein). The conservation and rea 



discussed in detail within the RMF framework 
symmetry in the Dirac negative energy spectrum 



proposed and investigated in RMF theory 



ization of pseudo-spin symmetry were 
With the same origin, the spin 
i.e. the single anti-nucleon spectrum) was 



As the Fock terms are missing and the one-pion exchange potential is not explicitly 



included in RMF, for the completeness of the theory, there have been attempts to inc 
the Fock terms in the ground-state energy of nuclear systems over the past two decades 



ude 

& 

20J. Recently, the RHF theory with density-dependent nucleon- meson couplings (DDRHF) 
finally succeeded in the quantitative description of the ground-state properties of many 
nuclear systems on the same level as RMF 2l|. Furthermore, it is found that the DDRHF 
theory can improve the descriptions of the nucleon effective mass and its isospin and energy 



dependences 



211 ]. as well as the shell evolution and closure with the inclusion of the one-pion 



exchange and p-tensor correlations [22], |23| . The pseudo-spin symmetry and its origin as well 



as the importance of the Fock terms have also been investigated before [241 . [25|. Although 
the pseudo-spin symmetry is still found to be a good approximation in RHF, its mechanism 
becomes rather complicated by the presence of the non-local potentials. 

In this paper, the Dirac negative energy spectrum or the single anti-nucleon spectrum 
in atomic nucleus such as 16 O will be investigated within the DDRHF theory in order to 
understand the relativistic symmetry with non-local potentials. The corresponding spin 
symmetry and its origin will be examined, in particular the role of the Fock terms. 
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II. THEORETICAL FRAMEWORK 



The starting point of the DDRHF theory is an effective Lagrangian density Jzf 21], which 
contains the degrees of freedom associated with the nucleon field (ip), two isoscalar meson 
fields (a and u), two isovector meson fields (tt and p) and the photon field (A). Then the 
effective Hamiltonian Jrf? is obtained with the general Legendre transformation. On the 
level of the mean field approximation, the energy functional $ is obtained by taking the 
expectation of the Hamiltonian , where both the Hartree (direct) and Fock (exchange) 
terms are kept. Finally, the Dirac equations, i.e. the equations of motion of nucleons, are 
obtained via the variation of the energy functional $ . 

For spherical nuclei, the nucleon Dirac spinor can be written as, 



fa(r) 



(1) 



where X|( r a) is the isospinor, ^ a a ma is the spherical harmonics spinor and ^j^ ma (v) = 
— (j ■ r^(j l ° ma (r) with l a = 2j a — l a . For the negative energy states, the lower component F(r) 
is dominant. The states are labelled by {nljm} with the relation 



n = n, for k > 0, 



n 



n + 1, for k < 0, 



(2) 



in analogy to Ref. 



The spin symmetry concerns the near degeneracy of the states 
(n, I, j — I ± 1/2). In the following equations, the sub-index will be omitted for simplicity. 
The radial Dirac equations are the coupled integro-differential ones due to the non-local 



Fock terms X and Y 
EG(r) = 

EF(r) -- 



18| 



d 
dr 



F(r) + [M + E s (r) + S (r)] G{r) + Y{r), 



+ 



d k> 
dr r 



G(r) - [M + S 5 (r) - S (r)] F(r) + X(r). 
Introducing the effective local potentials Xq, Xf,Yq and Yp by the definitions 

X(r) 



G<r)X #G(r) + fgf W . Xa (r)G(r) + X F (r)F(r) 



G 



G 2 + F 2 



(3a) 
(3b) 

(4a) 



G(r)Y(r) F{r)Y{r) 
K ' ' G 2 + F 2 K ' G 2 + F 2 



F(r) = Y G (r)G{r) + Y F (r)F(r) 



(4b) 
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the integro-differential equations Eq. fl3]) can be formally rewritten as equivalent differential 
ones, 

F(r) - [V+(r) - E] G{r) = 0, (5a) 



dr r 



d k 



+ X G (r 



dr r 

where V+ = VP + Y G , V. = V D + X F , and 



G(r) + [V-(r)-E]F(r) = 0, 



(5b) 



V° = M + E 5 + So, 



V D = E n - E,<? - M. 



(6) 



In the above expressions, E5 represents the scalar potential from the Hartree terms, E is the 
time component of the vector potential, which contains the contributions from the Hartree 
;he rearrangement terms induced by the density-dependence of the meson-nucleon 



terms and 

couplings [2l|, and X G , X F , Y G , Y F are the effective local potentials from the Fock terms. 
The equations Eq. (j3J) then can be solved self-consistently with the same numerical method 
as in RMF [271 . 



From the radial Dirac equation Eq. (jSJ), the Schrodinger-type equation for the dominant 
component F(r) can be obtained as, 

^ {F" + {V? + Vf) F' + [V CB + V^ P + Kf p] F} + V D F + V E F = EF, (7) 



V+ — E 



where Vcb 



and Vsop correspond to the centrifugal barrier (CB) and spin-orbit 



potential (SOP), respectively. In the above equation, the Hartree and Fock terms for V\, 
Vsop an d V read as 



Vi 



D 



vl 



D 

SOP 



V+-E' 

k vF' 



vf =x 



a 



V 



^=E 



rV + — E 

S5-M, 



yE 
v SOP 



V 



r 11/ 



V E =X F + 



E 
1 



V+-E' 
— X G — Y F 

Y. K 



— Y' F — X G Y F 



(8a) 
(8b) 
(8c) 



V+-E \ V + — E 

One may note that the denominator V + — E contains a state dependent potential Y G . How- 
ever, as the quantity Y G is around a few MeV and is negligible in comparison with V+ — E 
which is of the order of 1 GeV, the Eq. (ED) is accurate enough to estimate the Hartree 
and Fock contributions. Similar argument also holds for the time component of the vector 
potential Eq which contains the rearrangement term from Fock channels. 
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III. RESULTS AND DISCUSSION 



Solving the DDRHF equations Eq. ([5]) with the parameter set PKOl 21] in coordinate 
space self-consistently as in RMF j^], the neutron and proton single-particle energies can be 
obtained. We take the nucleus 16 O as an example to examine the negative energy spectrum 
and its spin symmetry. 
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FIG. 1: (color online) Single neutron spectrum in the Dirac sea for 16 calculated by DDRHF with 
PKOl. The dash-dot line represents the Hartree potential V D . For each pair of the spin doublets, 
the left levels are those with j = I — 1/2 and the right ones with j = I + 1/2. 



In Fig. [H all the bound single neutron states in the Dirac sea for 16 O are given. The dash- 
dot line represents the corresponding Hartree potential V D which is not strong enough for 
the 0s and Op orbits, the importance and contribution of the Fock terms is thus illustrated. 
For each pair of the spin doublets, the left levels are those with j = I — 1/2 and the right 
ones with j — I + 1/2. It can be clearly seen that the spin symmetry is well conserved in 
the Dirac sea. 

Taking p orbits with n — 0, 1, 2, 3 as examples, the Dirac wave functions of spin partners 
are shown in Fig. [2j The dominant components F(r) for the spin doublets are almost 
identical, whereas the small components G(r) show dramatic deviations from each other 
due to the node relation given in Eq. ([2]). The features of the spin partners for both the 
energies and wave functions are similar to those in RMF 15|. In the following, the origin and 
mechanisms of the spin symmetry will be investigated in comparison with those in RMF. 



The spin-orbit splittings AE\ S = E nl _ 



— Er in the negative energy spectrum of 
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O versus the average binding energies E av = {E n j_ + E n j, )/2 are given in Fig. [3j In 
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FIG. 2: (color online) Radial Dirac wave functions of the spin doublets p orbits in the negative 
energy spectrum of 16 O calculated by DDRHF with PKOl. Panels (a), (b), (c), and (d) are for 
Op, lp, 2p, and 3p spin doublets, respectively. 
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FIG. 3: (color online) Spin-orbit splitting AE\ S = E r —E, 



1+1/2 



J III 



1-1/2 



of 16 O versus the average binding energy E a 



III 



1+1/2 



+ E, 



ni 



1-1/2 



in the negative energy spectrum 
)/2 calculated by DDRHF with 



PKOl. The vertical dashed line shows the continuum limit. 



comparison with the RMF results (see Fig. 2 in Ref. 15j), the DDRHF results have the 
following characteristics: 1) the spin-orbit splittings are smaller; 2) the spin-orbit splittings 
fluctuate with E av , in contrast with the monotonous decreasing in the RMF case, when 
approaching the continuum limit; 3) in RMF the spin-down state (j = I — 1/2) is always 
lower than its spin- up partner (j = I + 1/2), while in DDRHF this occurs only for the p 
orbits and states near the continuum limit. 

In order to understand the origin of the spin symmetry in DDRHF and the relative 
positions of the spin-up state and its spin-down partner, the effective potentials V in Eq. 



6 



([7]) as well as the relations between the centrifugal barrier and the spin-orbit potential will 
be investigated. 




FIG. 4: (color online) Effective potentials in the negative energy spectrum of 16 calculated by 
DDRHF with PKOl, (a) for states np\/2 with n = 0, 1, 2, 3, (b) for states np^/2 with n = 0, 1, 2, 3, 
(c) for states 0pi/ 2 , 0d 3 / 2 , 0/ 5 / 2 , and 0<? 7/ / 2 , (d) for states 0p 3 / 2 , 0d 5 / 2 , 0/ 7 / 2 , and 0gg/ 2 . The Hartree 
part is labelled with dash-dotted lines. 

The effective potentials V for p, d, f, and g states in the negative energy spectrum of 16 
calculated by DDRHF with PKOl are shown in Fig. HJ together with the Hartree part V D 
(dash-dotted line). As seen in the Schrodinger-type equation Eq. ([7]), the effective potential 

V is composed of two parts, V D the Hartree potential from the direct terms, and V E the 
equivalent local potential from the exchange terms. The state dependence of the effective 
potential V comes from the contribution of the exchange terms. 

Corresponding to the nodes of the dominant component F(r), there exist fluctuations in 
the effective potentials V, which is brought in by the localization of non-local terms X and 

Y in Eq. ()4]). In addition, the contributions of Fock terms to the effective potentials tend to 
be slightly weaker when E av approaches the continuum limit, or for larger orbital angular 
momenta /. 

Comparing the left and right panels of Fig. HJ it is found that the effective potentials at 

r = are different between the spin partner states. This is due to the different asymptotic 

behaviors of the radial Dirac wave functions for spin doublets at r = 0, 

G(r) 

lim —. . oc r, for k > 0, 

F(r) (9) 



Fir. 

lim — . oc r, for k < 
r^o G{r) 



Within the RMF framework, it has been pointed out that the strong centrifugal bar- 
rier and weak spin-orbit potential lead to the pseudo-spin symmetry in the single nucleon 
spectrum 10J and the spin symmetry in the single anti- nucleon spectrum 15]. 





FIG. 5: (color online) Centrifugal barriers Vcb and spin-orbit potentials Vsop multiplied by the 
factor ^F 2 /(V+ — E) for the spin doublets (vOpi/2, v0p 3 / 2 ) (left panel) and (v^Pi^, v3p 3 / 2 ) (right 
panel) in the negative energy spectrum of 16 O. The insets show the Hartree contributions of the 
spin-orbit potentials. 

In Fig. [5] are shown the centrifugal barriers Vcb and the spin-orbital potentials Vsop 
multiplied by the factor ^fF 2 /(V + — E) for the spin doublets Op and 3p, and their integrals 
over r are respectively proportional to their contributions to the single-particle energy. It is 
clearly shown that the contribution of the centrifugal barriers Vcb is much larger than that 
of the spin-orbital potentials Vsop- Therefore, it can be concluded that similar reasons as 
in RMF lead to the spin symmetry in the negative energy spectrum in DDRHF, and the 
spin-orbit splitting is due to the different spin-orbit potentials Vsop of the spin doublets. 

In the insets of Fig. [5] are given the Hartree contributions to the spin-orbit potentials. 
It is found that the contributions from the Fock terms to Vsop are one order of magnitude 
larger than those from the Hartree terms. Therefore, the Fock terms must play important 
roles in the spin-orbit splitting of the spin doublets. 

From Eq. ([7]), the contributions to the single-particle energies E from different channels 
can be estimated quantitatively. For example, the CB contribution can be calculated by 



Jo 



F 2 dr 



V ( 



CB 



V+-E 



-F 2 dr. 



(10) 



In Table dare shown the contributions to the single-particle energies and spin-orbit splittings 



for the spin doublets Op and 3p. It is confirmed that the energy contributions from Vcb are 
much larger than those from Vsop and the contribution from the Fock terms Vggp is dominant 
in ^sop- 

TABLE I: The contributions from different channels (see Eq. (T7|)) to the single-particle energies E 
as well as the spin-orbit splittings AE for the spin doublets (vOpi/2, ^0p 3 / 2 ) and (v3pi/2, ^3p 3 / 2 ) 
in the negative energy spectrum of 16 0. The results are calculated by DDRHF with PKOl and all 
units are in MeV. 



state 


F" 


Vcb Vf 


^SOP 


V D 


Vf 


v SOP v 


E 


u0p 1/2 


-45.32 


-41.92 0.21 


-0.26 


-416.00 


5.35 


3.65 151.87 


-342.43 




-45.41 


-42.19 0.21 


0.53 


-415.68 


2.52 


8.18 149.49 


-342.37 


AE 


-0.09 


-0.27 0.00 


0.79 


0.32 


-2.83 


4.53 -2.38 


0.06 


v3p 1/2 


-211.40 


-38.44 0.01 


-0.12 


-611.42 


4.77 


0.90 55.82 


-799.91 


V3P3/2 


-211.39 


-38.42 0.01 


0.24 


-611.47 


3.80 


2.00 55.41 


-799.84 


AE 


0.01 


0.02 0.00 


0.36 


-0.05 


-0.97 


1.10 -0.41 


0.07 



In Table [H it is found that the contributions from Vg^ P , V E , Vggp and V E to the spin- 
orbit splitting are substantial. However, their contributions are counteracted by one another 
to preserve the spin symmetry. This kind of sophisticated cancellation implies that a weaker 
spin-orbit potential Vsop does not mean better conserved spin symmetry, as for Op and 3p 
orbits. 

To further confirm the role of the Fock terms, the contributions from the Hartree and 
Fock channels to the spin-orbit splittings in the negative energy spectrum of 16 versus the 
average energies of the spin doublets are shown in Fig. [6j where the Fock part includes the 
contributions from the terms Vf, V^p and V E and the rest is gathered into the Hartree 
part. It is found that the absolute contributions from both Hartree and Fock parts decrease 
monotonously with the average energy E av . The contributions from the Hartree terms have 
an energy dependence similar to those in RMF fisj ]. The contributions from Fock terms 
have an opposite tendency and cancel with the Hartree ones, thus leading to better spin 
symmetry. The competition between the Hartree and Fock terms will determine the sign of 
the spin-orbit splitting and this explains why the spin-orbit splittings in DDRHF fluctuate 
with E av in Fig. [3] instead of monotonously decreasing as in the RMF case. 
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FIG. 6: (color online) Hartree and Fock contributions to spin-orbit splitting in the negative energy 
spectrum of 16 versus the average energy of the spin doublets. The vertical dashed line shows 
the continuum limit. 

In order to get a deeper understanding of the cancellation and competition between the 
Hartree and Fock terms, we first separate the different meson contributions to single-particle 
energies. It is found that in both Hartree and Fock contributions, the isoscalar mesons, a 
and u, play dominant roles in the spin-orbit splitting, while the contributions from the p-, 
7r-mesons, and the rearrangement terms are negligible. Then, to make the mathematical 
structure simple and clear, one could replace the finite range Yukawa propagators with a 
pure delta function 5(ri — r 2 ) for these two heavy isoscalar mesons, but keeping their Dirac 
scalar and vector couplings. In this simple picture, it is found that the direct term of the 
cr-meson makes the spin-orbit splitting positive, whereas that of the w-meson makes the 
splitting negative. Since the attractive a field is slightly stronger than the repulsive u field 
in realistic nuclei, the net Hartree contribution to the splitting is positive, as shown in 
the upper part of Fig. [6l Meanwhile, it is also found analytically that the effect of the a 
exchange term is roughly one half as the effect of its direct term, but with an opposite sign. 
The effect of the u exchange term almost vanishes due to the cancellation between the time 
and space components. Therefore, the net Fock contribution to the splitting is negative and 
comparable to the Hartree contribution, as shown in the lower part of Fig. |6j All the above 
discussions for the case of 16 are also valid for heavier nuclei, e.g., 208 Pb. 

It is known that the presence of strong annihilation channels and various many-body 



effects cou 
couplings 



d_cause significant deviations from the G-parity values of the meson-antinucleon 
and a global fits to antiprotonic X-rays and radiochemical data indicates the 
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FIG. 7: (color online) Respectively same as Fig. Q] (a), Fig. [3] (b), and Fig. [6] (c), but with g^ff 

9iNN = 9iNN = Q-3diNN- 



need to use reduction factors [29|]. Therefore in order to exam whether the spin symmetry 
still persist in the face of these effects, the Dirac negative energy spectrum is calculated 
with the coupling constants g^N — 9inn = 9iNN — ®-3giNN, where i = a,ui,p. In Fig. 7, 
all the bound single neutron states in the Dirac sea thus obtained are shown, together with 
the spin-orbit splitting as well as its Hartree and Fock contributions. It is found that the 
spin symmetry is still well conserved even in the case where significant deviations from the 
G-parity values of the meson-antinucleon couplings due to the strong annihilation channels 
and various many-body effects occur. 



IV. SUMMARY 



In summary, the spin symmetry in the negative energy spectrum and its mechanism are 
investigated within the DDRHF theory by taking the nucleus 16 as an example. 

Similarly to RMF, the spin symmetry in the negative energy spectrum is found to be 
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a good approximation and the dominant components F(r) of the Dirac wave functions for 
the spin doublets are nearly identical, as the centrifugal barrier is much stronger than the 
spin-orbit potential. 

However, it is found that the Fock terms are dominant in the spin-orbit potential, which 
induce the state dependence of the effective potential and play essential roles in spin-orbit 
splitting. 

Classifying the contributions to the spin-orbit splitting into the Hartree and Fock parts, it 
is found that the Hartree part has an energy dependence similar to those in RMF 15| , while 
the Fock part has an opposite tendency and cancels with the Hartree part, thus leading to 
good spin symmetry. 

The competition between the Hartree and Fock terms determines the subtle spin-orbit 
splitting, which explains the fluctuation of the spin-orbit splittings in DDRHF instead of 
monotonously decreasing with E av as in the RMF case. 
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